Oscillatory Energy Exchange Between Waves Coupled by a Dynamic Artificial Crystal 
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We describe a general mechanism of controllable energy exchange between waves propagating in 
a dynamic artificial crystal. We show that if a spatial periodicity is temporarily imposed on the 
transmission properties of a wave-carrying medium whilst a wave is inside, this wave is coupled 
to a secondary counter-propagating wave and energy oscillates between the two. The oscillation 
frequency is determined by the width of the spectral band gap created by the periodicity and the 
frequency difference between the coupled waves. The effect is demonstrated with spin waves in a 
dynamic magnonic crystal. 



An artificial crystal is a wave transmission medium 
which features an "artificial lattice" created by a 
wavelength-scale modulation in its material properties. 
Such systems — which include, for example, optical pho- 
tonic crystals [3-0] and magnetic magnonic crystals 0- 
15| — belong to the class of so-called metamaterials: syn- 
thetic media with properties derived from an engineered 
mesoscopic structure. The wave transmission spectra of 
artificial crystals typically include band gaps. These gaps 
arise as a result of resonant wave-lattice interactions anal- 
ogous to the atomic-scale Brag g scattering phenomena 
observed in natural crystals [16| . Originally, the study of 
artificial crystals was directed exclusively toward time- 
invariant structures. However, considerable attention is 
now focussed on an emerging class of dynamic artificial 
crystals [H-H, H, 13 , 14 1 which have properties that can 
be modified whilst a wave propagates through them. 

A particularly interesting subset of dynamic artificial 
crystals have a lattice which can be turned "on" and 
"off" . Recently, it has been established that if such a 
crystal undergoes a rapid transition from "off" to "on" 
whilst a wave mode lying within the band gap is ex- 
cited inside it, this mode becomes coupled to a secondary 
counter-propagating wave, generally of a different fre- 
quency [14[. In this Letter, we demonstrate that the 
underlying physics of this coupling mechanism reveals 
a fundamental result of general wave dynamics. It is 
shown that the process is one of oscillatory inter-modal 
energy exchange. We explore the features of the effect 
experimentally using spin waves (magnetic excitations, 
the quanta of which are known as ma gnon s [T BI. Il7j ) in a 
dynamic magnonic crystal (DMC) [ullldl. Il5|. and show 
that results are consistent with a theoretical description 
of the coupling phenomenon which models the waves as 
interacting harmonic oscillators. 

Figure 1 (a) shows the experimental spin- wave system. 
The one-dimensional DMC comprises a planar metallic 
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FIG. 1: (Color online) (a) The dynamic magnonic crystal 
(DMC) comprises a metallic meander structure with 20 pe- 
riods of lattice constant a (10 shown) situated close to the 
surface of an Yttrium Iron Garnet spin-wave waveguide. Mi- 
crostrip input (left) and output (right) antennas are arranged 
on the surface of the film and a bias magnetic field Bq is 
applied along the spin-wave propagation direction, (b) Spin- 
wave transmission characteristics with the DMC in the "off" 
state (no current applied to the meander structure, dotted, 
blue) and in the "on" state (static current of 1 A in the me- 
ander structure, solid, red). In the latter case there is a band 
gap corresponding to the wavevector k a ~ n/a centered on 
Lo a = 2tt • 6500 MHz of -3 dB width approximately 30 MHz. 
(c) Calculated spin-wave dispersion relationship for the DMC 
system in the "off" (dotted, blue) and "on" (solid, red) states. 



meander structure (period a = 300 urn) situated close to 
the surface of an Yttrium Iron Garnet (YIG) thin-film 
spin- wave waveguide [H, [l5| . An antenna is provided at 
each end, and a static bias magnetic field Bo = 180 mT 
is applied along the spin-wave propagation direction. In 
the absence of a current in the meander structure, the 
transmission characteristics and spin- wave dispersion re- 
lationship uj(k) for the system are as shown by the dotted 
lines in Figs. 1(b) and 1(c), and the DMC is said to be 
"off" . If a current Iq (t) is applied to the meander struc- 
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FIG. 2: (Color online) (a) Inter-modal energy exchange in the 
DMC system. The diagonal lines in the upper panel represent 
the spin-wave dispersion curve. Solid dots (black) mark the 
wavevectors ±k a corresponding to the band-gap center fre- 
quency ui a ■ If the DMC is switched "on" at t = whilst an 
incident mode k 3 « k a , frequency lj s (black triangle) is inside, 
this mode is coupled to a secondary mode k r = k 3 — 2k a , fre- 
quency lu t (blue open circle). In the lower panel, the heights 
of the peaks centered on k s (black) and k r (blue) represent 
the power in the incident and secondary modes. Almost com- 
plete energy exchange occurs every T seconds for as long as 
the coupling persists, (b) Reversed (k r , blue open symbols) 
and transmitted (k s , black triangles) signal frequencies as 
a function of the applied spin-wave carrier frequency. Re- 
versed signal data is shown for two values of coupling dura- 
tion: r = 25 ns (circles) and r = 20 ns (squares). Transmitted 
signal datapoints for both values of r are represented by a sin- 
gle triangular symbol. Solid lines show the theory of Eq. (6). 
(c) Experimental reversed (P r , open blue circles) and trans- 
mitted (P s , black triangles) signal powers as a function of r 
for u) s = ui r — oj a - (d) Ratio of reversed to transmitted signal 
powers as a function of r; circular symbols (red) are points 
corresponding to the experimental data of (b), the solid curve 
(green) is the theoretical curve (see text). 



ture, its Oersted field (amplitude A£?o cx in) spatially 
modulates the waveguide's magnetic bias. Under these 
conditions, the DMC is "on" and has a band gap corre- 
sponding to the wavevector k a = n/a = 105 rad cm -1 
which appears as a discontinuity in the spin-wave dis- 
persion relationship centered on the resonance fr equency 
uj{k a ) = oj a = 2tt-6500 MHz (solid line, Fig. 1(c)) U^M- 
The band-gap width is determined by ABq (~ 1 mT in 
our experiments) and thus the current amplitude. 

If a spin wave with wavevector k s w k a = ir/a is in- 
cident on the crystal whilst it is in the "on" state, it 
is simply reflected, creating a "dip" in the experimen- 
tal transmission characteristic (Fig. 1(b)). However, if 



the crystal is switched from "off" to "on" whilst such a 
wave is inside, a geometrically-induced coupling between 
the incident mode and secondary modes with wavevec- 
tors k r — k s ± 2k a is brought about. This coupling 
is pronounced when the eigenfrequencies of both waves, 
lo s = tj(k s ) and u r = U)(k r ), are close to each other, 
i.e. for k s k a and thus counter-propagating waves 
k r = k s — 2k a w — k a . The dynamics of the complex 
amplitudes c s and c r of the two interacting waves can be 
modeled by a pair of coupled equations [19( 



dc s 
~dt 
dc r 

~dt 



iuj r C r + ifl a C s . 



(la) 
(lb) 



where Q a is the coupling constant (with units of angular 
frequency), numerically equal to the half- width of the 
crystal band gap. Oscillatory solutions to (la) and (lb) 
have frequencies 



;i )2 = 57 ± v^a + ^ 2 = 57 ± Q 



(2) 



where 57 = (ui s + uj r )/2 is the average eigenfrequency 
of the waves, and Soj = (lu s — cu r )/2 is their frequency 
mismatch. Equation (|2|) describes the dispersion relation 
for a static magnonic crystal in the "on" state with a band 
gap of width 2f2 a centered on the resonance frequency 

LU S = UJ r = UU a . 

In a dynamic experiment, at the moment when the 
crystal is switched "on" (t — 0) the amplitude of the 
incident wave is non-zero (c s = 1), whilst the coupled 
secondary wave is absent (cy = 0). Solving Eqs. ([T]) 
with these initial conditions, one can find how the power 



p — 



of each wave changes with time: 



8lo 2 

PJt) = cos 2 (fit) + — sin 2 (Ot) , (3a) 
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(3b) 



The equations above describe periodic energy ex- 
change at frequency f2 (see Eq. ([2])) between confined, 
geometrically-coupled incident (k s ) and secondary (k r ) 
modes. The maximum power fraction transferred to the 
k r mode R max = ^ 2 /^ 2 = ^l 2 a /{^l 2 a + Suj 2 ) is of order 
unity for waves within the band gap (\Su>\ < Q a ) and is 
achieved after a characteristic energy exchange time 



T = — = = 

2^ 2yjni + 5uj 2 



(4) 



or odd integer multiplies of thereof (see Fig. 2(a)). 

If the dynamic crystal is switched "off" at t = r, elim- 
inating the coupling, the waves, which have powers given 
by the instantaneous solutions of Eqs. ([3]) travel in op- 
posite directions. Experimentally, the oscillatory energy 
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exchange process may be studied by recovering the ratio 
of the powers of these propagating waves, 



«(r) 



P r (r) 



n 2 a tan 2 (Or) 



P s (t) W + Soj 2 tan 2 (fir) ' 



(5) 



which is independent (up to a constant factor) of their 
damping. Note that the coupling effects frequency 
conversion — transfer of energy from the incident to 
the secondary mode — through an entirely linear mech- 
anism 14]. Since this process is efficient only over a 
narrow range of wavevectors k s s=s k a , the dispersion re- 
lation can be approximated by a first order Taylor ex- 
pansion, lo s = io(k s ) = uj a + v g (k s — k a ), uj r — u)(kr) — 
oj a — Vg(k r + k a ), where v g = dio(k)/dk is the group veloc- 
ity [2 Of at k = k a , and we assume a symmetric dispersion 
law ui(k) = u)(—h): 



UJ r — 2u) a 



(6) 



i.e. the detuning of the incident and secondary mode 
frequencies is inverted with respect to uj a . 

Experimental investigations into the geometrical mode 
coupling mechanism were performed using a current 
pulse Io (t) of duration t to briefly switch the DMC "on" 
whilst an incident spin-wave signal packet of center fre- 
quency uj s and duration 200 ns excited by a microwave 
drive signal (An, Fig. 1(a)) was contained within it 21 1. 
When the DMC is switched "on" (t = 0), transfer of 
energy from the incident k s mode into the secondary k r 
mode initiates (see Fig. 2(a)), both modes being confined 
to the crystal. Following the deactivation of the coupling 
at t = t, the k s mode propagates to the output antenna, 
producing a transmitted signal of power A s oc P s (r), 
whilst the k r mode travels to the input antenna, orig- 
inating a "reversed" signal of power A r oc P r (T~) (see 
Fig. 1(a)). (Note that, in the absence of damping and 
with perfect antenna efficiency A s>r — P s ^(t)). 

Measured frequencies of the transmitted (k s ) and re- 
versed [k r ) signals are plotted in Fig. 2(b) for two values 
of the coupling duration r. The frequency of the trans- 
mitted signal is not changed by interaction with the crys- 
tal. The reversed frequency, as predicted by Eq. ©, is 
determined solely by the geometry of the DMC and is 
independent of r. The data are a near-perfect fit to the 
theory (solid lines on the same axes), providing strong 
evidence for its validity. The antiphase variation of P s 
and P r with the coupling time r, predicted by Eqs. (J]), 
is confirmed by the experimental data of Fig. 2(c) for the 
resonance case to s = u> r = oj a . Figure 2(d) shows the ratio 
of the reversed to transmitted signals of Fig. 2(c) (circular 
symbols) and the corresponding theory (k — P r /P s , solid 
line). As can be seen, there is good agreement between 
the two trends. For the theoretical curve of Fig. 2(d) we 
assume fl a = 27T-5.53 MHz, a small detuning Slo = 0.03%, 
and an effective coupling duration r e s = r — 2Ar where 
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FIG. 3: (a) Reversed (A r oc P r (j), Fig. 1(a)) signal spectra 
for varying r (input spin-wave signal duration 150 ns, carrier 
frequency ui s = u) a = 2ir ■ 6500 MHz). Red is indicative of the 
highest signal intensities, dark blue the lowest, each spectrum 
is individually normalized, (b) Theoretical calculation of the 
data of (a) using experimentally derived parameters. Panels 
(c) and (d) show sections of the experimental data of (a) (red 
circular symbols) and (b) (blue curves) at r = 50 ns and 
r = 105 ns respectively (vertical dashed lines, (a) and (b)). 



2 At = 12.5 ns to account for the effect of the finite rise 
and fall times of the DMC current pulse. The plot is also 
scaled to match the experimental and theoretical max- 
imum ratios (absolute values being dependent on losses 
and antenna efficiencies). The reduction in the signal 
amplitudes of Fig. 2(c) with increasing r occurs as a re- 
sult of damping. The non-zero minima of the reversed 
curve results from the finite spectral width of the input 
signal (which accounts too for the difference in heights 
of the two peaks in Fig. 2(d), the total power of the re- 
versed signal is the integral of Eq. (3b) over all frequency 
components of the input signal) and weak spurious reflec- 
tion of spin waves from the metal meander structure. In 
combination with damping, spurious reflections are also 
responsible for the small offset between the maxima and 
minima of the transmitted and reversed curves. 

Note that the value of the coupling constant f2 a = 
2tt ■ 5.53 MHz derived from the dynamic experiments re- 
flects the width of the static DMC band gap at approx- 
imately —20 dB (Fig. 1(b). This is consistent with the 
fact that the coupling mechanism requires strong rejec- 
tion of waves from the crystal. 

For a fixed incident mode power, Eqs. (3) predict that 
the frequency dependence of the secondary mode power 
varies with the coupling duration t. In order to observe 
this dependence experimentally, the power spectrum of 
the reversed signal (A r (t), Fig. 1(a)) produced from an 
input spin-wave packet of duration 150 ns and carrier fre- 
quency u> s = uj a = 2ir ■ 6500 MHz was recovered for vary- 
ing t. Results are plotted in Fig. 3(a). Figure 3(b) shows 
the corresponding theory. As can be seen, the theory ac- 
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FIG. 4: Theoretical spatio-temporal maps of interacting 
waves in DMC system for different coupling durations r at 
oj s = oj r = ui a . Each panel shows the space- (horizontal) time 
(vertical) dependence of the wave power. Red and blue colors 
correspond to the incident (k„) and secondary (k r ) waves re- 
spectively and the intensity of the color represents the wave 
power. Dashed horizontal lines indicate the time interval r 
during which the DMC is "on" . (a) r = T. (b) r = 2T. (c) 
r = 3.5T (50:50 beamsplitter), (d) r -> oo. 



curately predicts the key features of the data, namely the 
high signal intensity at ui r = u a which extends from the 
shortest measured coupling duration (15 ns) to around 
t = 90 ns, and the branching of the spectra at approxi- 
mately t — 90, 180 ns. Figures 3(c) and 3(d) are sections 
through the data of Figs. 3(a) (red circular symbols) and 
3(b) (blue curves) at r = 50, 105 ns respectively. The de- 
viation between theory and experiment at long coupling 
times in Figs. 3(a) and 3(b) occurs as a result of leakage 
of energy from the DMC during the coupling interval; 
an inevitable feature of any experiment using an artifi- 
cial crystal of finite length. The different peak heights in 
Fig. 3(d) are caused by a small detuning (0.03 %) of the 
carrier signal frequency from the resonance frequency u) a 
in combination with a weak frequency dependence of the 
characteristics of the microstrip antennas. 

To further clarify the features of the mode coupling 
mechanism, Fig. 4 shows theoretical spatio-temporal 
maps of the interacting waves in the artificial crystal sys- 
tem for different coupling durations r. In each panel, 
the power of the incident (k s , red) and secondary (fc r , 
blue) modes are plotted as a function of time (vertical) 
and position (horizontal) and the dashed horizontal lines 
indicate the time interval r. An ideal, lossless system 
with uj s = ijj T — uj a is assumed. In the limit r — > oo the 
incident wave energy is effectively "stopped" ; it remains 
confined to the crystal and repeatedly oscillates between 
the k s and k r modes (panel (d)). 

In summary, we have described a mode coupling phe- 
nomenon in a dynamic artificial crystal. The mechanism, 
which is a principle of application to waves of any na- 
ture, involves oscillatory energy exchange between two 
counter-propagating wave modes, generally of different 
frequencies. The repeated spontaneous reversal of energy 



flow direction which characterizes the effect distinguishes 
it from all other known methods of inter-modal e nerg y 
transfer in wave systems (see for example [HI, [13, ) . 
Experimental measurements performed using a dynamic 
magnonic crystal are in agreement with a theoretical de- 
scription of the mechanism which models the coupled 
waves as interacting harmonic oscillators. The relation- 
ship between the frequencies of the coupled modes is de- 
termined by the crystal geometry, and the rate of energy 
exchange between them by the width of the crystal band 
gap and their frequency mismatch. As well as represent- 
ing a fundamental development in the understanding of 
wave dynamics in metamaterial systems, the effect po- 
tentially finds widespread technological application in the 
context of signal processing. 
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